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A technique is proposed for creating nonground-state Bose-Einstein condensates in a trapping 
potential by means of the temporal modulation of atomic interactions. Applying a time-dependent 
spatially homogeneous magnetic field modifies the atomic scattering length. A modulation of the 
scattering length excites the condensate, which, under special conditions, can be transferred to an 
excited nonlinear coherent mode. It is shown that a phase-transition-like behavior occurs in the 
time-averaged population imbalance between the ground and excited states. The application of the 
technique is analyzed and it is shown that the considered effect can be realized for experimentally 
available condensates. 

PACS numbers: 03.75.Kk,03.75.Lm,03.75.Nt 
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I. INTRODUCTION 



Bose-Einstein condensation (BEC) is an effect of a very 
broad interest, touching on a variety of physical topics 
and creating, sometimes unexpected, links among differ- 
ent physical subjects. Many exciting possibilities have 
been investigated in recent years and are discussed in 
books [l| and review articles [1, S 3 ■ The Feshbach reso- 
nance technique, which enables a variation of the scatter- 
ing length via a magnetic field, is one of the most promis- 
ing tools for manipulating the properties of quantum de- 
generate gases. This technique, e.g., has made it possible 
to tune an unstable system into a stable one 0,0], to form 
molecular condensates [3], and to investigate the nonlin- 
ear dynamics of a BEC with a time-dependent scattering 
length gg. 

In the present letter, we show that the temporal mod- 
ulation of the scattering length can also be used for gen- 
erating nonground-state condensates of trapped atoms. 
Such states are described by nonlinear topological coher- 
ent modes and can be excited by a resonant modulation 
of the trapping potential El • We advance an alter- 
native way for exciting the coherent modes of a trapped 
BEC by including an oscillatory component in the scat- 
tering length. The main idea is to superimpose onto the 
BEC a uniform magnetic field with a small amplitude 
time variation. Due to the Feshbach resonance effect, 
such an oscillatory field creates an external perturbation 
in the system, coherently transferring atoms from the 
ground to a chosen excited coherent state. The feasibility 
of the experimental implementation of this phenomenon 



for available atomic systems is demonstrated. 

The coherent states of a trapped BEC are described by 
the solutions to the Gross-Pitaevskii equation (GPE). To 
transfer the BEC from the ground to a nonground state, 
it is necessary to apply a time-dependent perturbation, 
at a frequency close to the considered transition. As a 
result [m the resonantly excited condensate becomes 
an effective two-level system. The external fields consid- 
ered in the previous works 0, [13, [O, were formed by 
spatially inhomogeneous alternating trapping potentials. 
Now, we consider a very different situation represented by 
a spatially homogeneous time-oscillating magnetic field, 
which can be easily implemented with present experimen- 
tal techniques. 



II. MODULATION OF SCATTERING LENGTH 

The GPE, describing a zero-temperature weakly inter- 
acting Bose gas, is given by 
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where Utrap (r) is the trapping potential, the interaction 
strength Ag — 4:n{N — l)fi'^as/mQ, the s-wave scatter- 
ing length. Too the atomic mass and N is the number of 
condensed atoms. In the presence of a spatially uniform 
magnetic field, near a Feshbach resonance is given by 
the well known relation 
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where a„r is a non-resonant scattering length, Bres is 
the value of the magnetic field where the resonance in 
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occurs, and A is the corresponding resonance width. Let 
us consider the time-dependent magnetic field 



B{t) = Bo + bcos{Lut), 



(3) 



with l&l ^ \Bo — Bres\- In such a case, Eq.(l2|) can be 
expanded to first order as 



as{t) ~ flQ + acoa{ujt), 



where 



Bo — Br 
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The scattering length then possesses an oscillatory com- 
ponent around the average value. 

Combining Eq.dU and Eq.([T]), one gets the GPE with 
the additional oscillatory term V ~ V {r,t). With the 
notation H — Ho + V, one has 



i7$ = iHd<^>/dt, 



(6) 



Ho = -^V^ + Utrap{r) + Aom\ (7) 



V = Acos{Lut) 1$ 
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Solving Eq. ([6|) , we keep in mind that the frequency uj is 
chosen to be close to the transition frequency between the 
ground and an excited mode. We starting considering as 
the total wavefunction a linear combination of a complete 
set of modes as follows 



$ (r, t) 
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where 4'j(r) are stationary solutions for the equation 
Ho4>j = Ej^ji with eigenenergies Ej. Then, was proved 
before 0, the only relevant terms that survive are two 
modes connected by the modulating perturbation ([8]) . In 
this case, the total wavefunction ^ can be represented, 
in a good approximation, by 

$ (r, t) = co(t)0o(r)e-^^°*/'^ -f Cp(t)0p(r)e-'^-'/^ (10) 

where the label refers to ground state and p to an 
excited state. We investigate the time evolution of the 
BEC, with the initial condition, where all atoms are in 
the ground state, i.e., Co(0) = 1 and Cp(0) — 0. Our 
aim is to study the transfer between the fractional mode 
populations of the ground and the excited states. Using 



Eq. pO]) in Eq.®, one obtains a set of differential equa- 
tions for the coefficients co{t) and Cp(i), 
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where the integral Ij^k.i is defined as 



(f>*\(t)k\^(f)idr. 



(lib) 



(12) 



In deriving the latter equations, two assumptions, whose 
mathematical basis has been described in detail in 
Refs. 0, Hi Hill are made. First, the time variation 
of Co(t) and Cp(t) are to be much slower than the the 
exponential oscillations with the transition frequency 
'-^po = {Ep — Eo)/h. This condition is fulfilled, when 
the amplitudes AoIj^k,i and Alj^k.i are smaller than hujpo. 
The second is the resonance condition, when the external 
alternating field connects only the two chosen nonlinear 
states. Another point concerns damping due to colli- 
sions between particles in the desired modes or collisions 
with the thermal cloud. Although the oscillation time for 
populations takes tens of trap periods, this time is much 
smaller than the lifetime of a typical BEC or a vortex 
state [H, [l2|- So, we expect that damping occurs but 
not as a dominate process. Thus, we have left out the 
damping effect for this model. 

Another important aspect is that the total number of 
atoms do not vary on time, but the number in each state 
does. This variation is taking into account in Eqs. pT|) 
since these equations depend on the population of each 
(9) state, represented by |coP and |cpp. However, the modes 
in equation ^ are stationary solutions of Equation 
([7]) when all atoms are in state j. Thus, if there is a 
variation in the atom number of some state, there is a 
variation in the wave function that represents this state. 
So, the total wave function should be written in the form 



cl>(r,i) = ^d,(t)0;(r,i). 



(13) 



where the number dependence is inserted in the time de- 
pendence. In this way, the population of a state j would 
be given by |dj(f)p and not by |cj (t)p, since the expan- 
sions (III) and ([T3|) are different. However, in the case of 
our study, the system is in a weak-coupling regime, i.e., 
g is small, so the variation of the wave function can be 
neglected and the population of a state j can be given by 
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III. APPLICATION TO A CYLINDRICALLY 
SYMMETRIC TRAP 

We consider a cylindrically symmetric harmonic trap 



TT / 2 2 I 2 2 



(14) 



and use the o ptim ized perturbation theory, as discussed 
in Ref. [3. ITol.TTl| . for finding the modes 0o and (j)p. It is 
convenient to define the dimensionless variables 



^(x) =/3/2^(r), = 



tklJr 
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go = An{N - 1 



A = 



S=t^^, (15) 



where Ir = 
Kutta method 
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rnoUJr- Using the fourth-order Runge- 
we calculate the time evolution of 
the coefficients co{t) and Cp{t) for different values of the 
detuning 6 and external-field amplitude. The functions 
co(t) and Cp{t) define the fractional mode populations. 

For an excited mode, we take the radial dipole state 
{100}, which is the lowest energy state above the ground 
state {000}. Here {nmk} implies the notation for the 
three quantum numbers, where n, m and A; refer to radial, 
azimuthal and axial mode numbers respectively. FiglU 
where A = 0.2 and go = 70, shows the time evolution of 
the mode populations hq and Up for different values of the 
detuning S and a/ao, which is given by Eq. ([5]). The cho- 
sen parameters correspond to typical experimental setups 
and are easily controlled in a laboratory. The solutions 
demonstrate different behaviors of the state populations. 
For a/ao = 0.7 and (5 = 0, Fig[TJa), the populations 
display small oscillation amplitudes, with a considerably 
larger population in the ground state. Increasing the 
detuning to 6 = 0.04 results in the behavior shown in 
FiglUb). Although on average atoms stay longer in the 
ground state, for some intervals of time Hp is larger than 
no- Changing the amplitude to a/ao — 0.72 and main- 
taining S — 0.04, as in Fig[TJc), yields a very different 
temporal behavior. Atoms now stay longer in the ex- 
cited state rather than in the ground state. The shape of 
the functions shows the inherent nonlinearity of the sys- 
tem. If the amplitude is increased further to a/ao = li 
with 6 — 0.04, as in Fig. [Ijd), the system shows a full 
population inversion. For certain times, when the mode 
population fully migrates from the ground {000} to the 
excited {100} state, it is possible to have a pure conden- 
sate in the coherent topological excited mode. 

A convenient way to quantify the population behav- 
ior is through the introduction of an order parameter ry, 
defined as the difference between the time-averaged pop- 
ulations for both states [13, [ll| , 




r] = no- Up. 



(16) 



t (units of co"^ ) 



FIG. 1: (Color online) Populations of the ground state no 
(solid line) and excited state {100} Up (dashed line) as a func- 
tion of time for A — 0.2 and go ~ 70 with (a) a/ao — 0.7 and 
(5 = 0; (b) a/ao = 0.7 and S = 0.04; (c) a/ag = 0.72 and 
S = 0.04; (d) a/ao = 1 and 5 = 0.04 



Here, the average of each population is performed over 
the full cycle of an oscillation. 

The above order parameter rj displays a nontrivial be- 
havior when the ratio a/ao is modified. For different de- 
tunings, the variation of 77 as a function of the ratio a/ao 
is presented in Fig[2l The variation of 77 vs. a/ao can 
be smooth, when S < 0.03, or can show sudden changes, 
when S > 0.03. Smaller values of a/ao keep atoms prefer- 
ably in the ground state. For some critical value of a/ao, 
rj becomes negative, which means that BEC stays longer 
with a larger population in the excited state than in the 
ground state. This situation is ideal for detecting the 
formation of topological modes. Step like behavior ob- 
served on 77 is not surprising for a nonlinear system. For 
example, the classical driven anharmonic oscillator [iTj 
exhibits a shift on its resonance curve, which we also ob- 
serve in our system [18[. This shift for a large driven 
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FIG. 2: (Color online) Order parameter 77 as a function of the 
ratio a/ao, for different detunings showing a phase-transition 
like behavior. 



strength can result in a bistable behavior. In our system 
that appears as a step like response with ?/. This bista- 
bility appears only for positive detunings, as observed 
in our calculation. Negative detunings follows a smooth 
curve without step occurrence. 

An important question is the feasibility of the exper- 
imental creation of such coherent modes. The main pa- 
rameter here is a/ao. It shows us that this value is 
strongly dependent on the Feshbach resonance width A, 
characteristic of each type of atoms. For systems with 
small A, the required value of a/ao for obtaining the 
transition in ry will occur only for Bq w Bres- In this 
case, the necessary condition b <C Bres — Bq can only 
be fulfilled for very small values of &, which creates an 
extra difficulty with the present techniques of magnetic 
field control [19, 20J. As an experimentally realistic ex- 
ample, let us consider the case of 6 = 0.1{Bres — Bq) and 
a/ao = 0.8, which corresponds to the atomic parameters 
listed in Table [H Setting go = 70, = 27r x 120 Hz, 
we obtain a value of b of 0.9 G for ^^Rb, and of 0.02 G 
for ®^Rb, which would be difficult to control. On the 
other hand, 10 G, for ^Li, and 4.55 G, for '^^K, are val- 
ues that can be realized with present technical capabili- 
ties til, [23. 

Let us consider now a condensate containing 10^ ^Li 
atoms in a trap with radial frequency tOr — x 120 Hz, 
go = 70, and A = 0.2. With these conditions, the tran- 
sition frequency is ujpo = wioo,o ~ ^tt x 209.6 Hz. Also, 
together with the information from Table |T1 Eqs. (O and 
P^ . we obtain the bias magnetic field Bo — 632.5 G 
and the amplitude b — 11.68 a/ao. Then, we observe 
the critical values for b ranging from 4 G for a detuning 
uj — ujpo = 27r X 21 Hz to 10.5 G, for cu — uipo — 27r x 6.3 Hz. 
Such oscillating amplitudes correspond to less than 1% 
of the total bias field Bq. 



Atom 


Bres {G) 


A(G) 




Bo(G) 


6(G) 


ao 


iV(xlO^) 


«5Rb 


155.0 


10.7 


-443 


164.4 


-0.9 


63 


0.2 


»^Rb 


1007.34 


0.17 


100 


1007.53 


0.02 


11 


0.9 


^Li 


735 


-113 


-27.5 


636 


10 


3.9 


9.3 


39^ 


403.4 


-52 


-23 


357.9 


4.55 


3.3 


4.7 



TABLE I: Amplitudes Bo and 6 of the magnetic field for 
four difi^erent species of atoms . We set go = 70, A = 0.2, 
a/ao = 0.8 and b — 0.1(-Bres — Bq). Scattering length is 
expressed in units of the Bohr radius. Data for *^Rb were 
obtained in Ref. for ^'^Rb, Ref. (22]; for "^Li, Ref. [3; for 
^^K, Ref. M. 



A final point to be addressed concerns losses intro- 
duced by collisions, specially near a Feshbach resonance. 
With an off-resonant magnetic field, the dominant loss 
mechanism is a three body collision, whereas close to 
the resonance, the molecular formation dominates the 
atom loss mechanism [2^. Although the fields consid- 
ered in Table U and in the ''Li example above are within 
the resonance linewidth, they are far enough to the res- 
onance, then we consider three body collision as the 
main loss mechanism. In the case of ^Li, that the res- 
onance linewidth is large, as we are in the border of 
this linewidth, the loss rate is very small, as shown in 
Ref. [ll]. So, we can neglect the atom loss, specially be- 
cause we apply the magnetic field in a short period of 
time. 



IV. DISCUSSION 

In conclusion, we have shown that using a magnetic 
modulation field, applied to a trapped Bose-condensed 
gas, it is possible to transfer the atomic population 
from the ground state to an excited state, producing a 
nonground-state condensate. The time averaged popu- 
lation imbalance between the ground and excited states 
represents an order parameter, which demonstrates an 
interesting behavior as a function of the modulation am- 
plitude. Depending on the detuning, the behavior of 77 
can be either smooth or rather abrupt. This is the conse- 
quence of the strong nonlinearity of the interactions. For 
some range of detunings, rj becomes negative above a crit- 
ical value of the modulation amplitude. This occurs be- 
cause of the population inversion realized during the pro- 
cess of the mode excitation. Larger detunings and out-of- 
resonance modulations keep the population in the ground 
state and no population inversion is observed. Numerical 
calculations, accomplished for ^Li atoms, show that the 
values for the amplitude and modulation of the bias field 
are within realistic experimental conditions. 

The justification of the approach for generating non- 
linear coherent modes by the resonant modulation of an 
alternating field has been explained in detail in our pre- 
vious publications [H, Uil, Hill HI- It would be unrea- 
sonable to repeat here all this justification in full. How- 
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ever, for the convenience of interested readers, we recall 
in brief some important points. 

When a system is subject to the action of a time- 
dependent external field V{t), alternating with a fre- 
quency w, then there exist two principally different situa- 
tions, depen ding on the temporal alternation being either 



26[. Respectively, there occur two different 



slow or fast 
physical cases. 

If the temporal variation is slow, such that the alterna- 
tion frequency uj is much smaller than the characteristic 
system energy E, then one says that the perturbation is 
adiabatic [20|. For a quantum- mechanical system, the 
energy E is an eigenvalue of the system Hamiltonian H . 
When the latter varies slowly in time, so that the varia- 
tion is adiabatic, then the adiabatic picture is in order, 
with the eigenvalues E{t) and related eigenfunctions ip{t) 
slowly varying in time. For such an adiabatic variation, 
the notion of nonlinear coherent modes is not of much 
use. 

A principally different situation occurs, when the tem- 
poral variation is fast, so that uj is in resonance with 
one of the transition frequencies LUij = {Ei — Ej)/h. 
It is exactly this case that is considered in the present 
paper. Then the solutions of the stationary equation 
Hoipj = Ejipj define the nonlinear coherent modes ipj. 
Constructing a closed linear envelope, spanning the total 
set {<fj}, one gets the Hilbert space H = Span{ipj}, with 
the typical property [13] that any function from H can be 
expanded over the set {(fij}- The latter set forms a total 
basis [ni nil J over which the solution of Eq. (6) can be 
expanded as in Eq. (9). Looking for the solution of Eq. 
(6) in the form of expansion (9), we employ the method 
of the parameter variation. Separating in expansion (9) 
the time-dependent coefficient function into two factors, 
the fastly oscillating exponential exp(—iEjt/h), and the 
slowly var ying envelope Cj(t), we can use the averaging 
techniques [1^ and the scale separation approach [lOiUll . 
A similar technique in optics is called the slowly-varying 
amplitude approximation [32]. Substituting expansion 
(9) into Eq. (6) and involving the averaging techniques 
[29-31], we obtain Eqs. (11) for c„(t). This procedure is 
the same as has been done in Refs. [10,11]. 

In this way, one should not confuse the adiabatically 
slow variation of an external field, when the mode profiles 
and energies are certainly changing in time, and the fast 
resonant field oscillation, when the fractional mode pop- 
ulations vary between the stationary coherent modes. In 
the latter case, the solution to Eq. (6) can be represented 
in form (9) and the averaging techniques are applicable, 
yielding Eqs. (11) for the coefficient functions. The sit- 
uation here is equivalent to the resonant excitation of 
an atom, as is discussed in Refs. [24,25]. In the same 
way as for an atom, the external resonant field should 
not be too strong in order not to disturb the energy-level 
classification. For this, it is sufficient to take the modu- 
lation amplitude b in Eq. (3) appropriately small, which 
is always possible. 

It is important to stress that the averaging technique, 



employed for deriving the equations for the fractional 
mode coefficients c„(i), is a well justified method based 
on rigorous mathematical theorems [29]. Therefore, this 
method does not require some other justifications. 

Moreover, the results, obtained by means of the av- 
eraging techniques, when treating the generation of the 
nonlinear coherent modes by a resonant alternating field, 
have been thoroughly compared with the results of the 
direct numerical simulation of the Gross-Pitaevskii equa- 
tion in Refs. [33,34]. Both ways of treating the problem 
have been found to be in a very good quantitative agree- 
ment. 

In a real experiment, of course, a modulation, even be- 
ing perfectly resonant, might, nevertheless, lead to heat- 
ing, when the neighboring nonresonant levels become es- 
sentially involved in the process. This means that such 
a resonant mode generation can be effective only during 
a finite time. This problem has been analyzed in detail 
earlier [11,35], where it has been found that, taking the 
modulation amplitude sufficiently small, it is feasible to 
avoid heating during the lifetime of atoms inside a trap. 
For instance, the heating time was estimated [11,35] to 
be of order 10 - 100 s. 

The nonground-state nonlinear coherent modes are of- 
ten termed topological, since the wavefunctions of dif- 
ferent modes possess drastically different spatial shapes, 
with a different number and location of zeros. As an 
example, we can recall the wavefunctions of the first ex- 
cited modes for atoms in a harmonic trap, found in Refs. 
[11,24]. These functions ipnmkixrjfjXz), represented in 
dimensionless units of Eq. (14), are labelled by the quan- 
tum numbers n, m, and k and are the functions of the 
dimensionless cylindrical variables Xr,ip, and Xz- The 
ground-state wavefunction is 

/ 2 \ 1/4 

The radial dipole mode has the form 

V^ioo = (%^)'^' {uiooxl - 1) e-(-«-^+-«-^)/2 . 
The vortex mode is also a possible stationary solution 



And the axial dipole mode is 



The quantities u,mifc and Vnmk here are defined by the 
optimization conditions and depend on all system pa- 
rameters (see calculational details in Refs. [11,24]). As 
is evident, the spatial shapes of atomic clouds, described 
by ]V^nmfc]^7 are principally different for different quan- 
tum numbers. 
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If the scattering length is modulated according to Eq. 
(4), then the resulting Eq. (6) acquires the additional 
term (8) playing the role of a modulating field. Therefore 
the mechanism of generating nonground-state modes is 
the same for both the cases, whether the additional term 
is caused by the trap modulation or by the scattering- 
length modulation. In any case, the general structure of 
Eqs. (6) to (8) is similar for both these setups. 

The main difference between the methods of generat- 
ing nonlinear coherent modes by the trap modulation or 
by the scattering-length modulation is as follows. In the 
method of the trap modulation, by choosing the appro- 
priate spatial dependence of the modulating field V{r, t), 
it is possible to generate any type of modes. While, in the 
method of the scattering-length modulation, only those 
modes can be excited, for which the integrals in Eq. (12) 
are nonzero. For example, in the case of cylindrically 
symmetric harmonic trap with potential (13), the inte- 
grals Io,p,p, Iofi,p, Ip,o,o, and Ip,p,o, in which p impUes a 



set {n, m, k}, vanish for the vortex mode ipp = "0010 and 
for the axial dipole mode ipp = V'ooi- Hence, in these 
cases, the vortex and the axial dipole modes cannot be 
excited by modulating the scattering length. However, 
the radial dipole mode tpioo can be generated by this 
method, since the related integrals (12) are nonzero. This 
is why, we have considered here exactly this case. A way 
out of the problem, which would allow for the generation 
of other modes, including the vortex mode, could be by 
employing a nonsymmetric trapping potential. 
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